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Abstract 

A one-loop effective potential of the Higgs field on the Schwarzschild background is derived in the frame- 
work of a toy model: a SO(N) scalar multiplet interacting with the gauge fields, the SO(N) gauge symmetry 
being broken by the Higgs mechanism. As expected, the potential depends on the space point and results 
in a mass shift of all massive particles near a black hole. It is shown that the obtained potential depends 
on the space point through the metric component goo in the adapted coordinates and has the same form for 
an arbitrary static, spherically symmetric background. Some properties of this potential are investigated. 
In particular, if the conformal symmetry holds valid for massless particles on the given background, there 
exist only two possible scenarios depending on a sign of an arbitrary constant arising from the regularization 
procedure: the masses of all massive particles grow infinitely when they approach the black hole horizon, or 
the gauge symmetry is restored at a finite distance from the horizon and all particles become massless. If the 
conformal symmetry is spoiled, an additional term in the effective potential appears and the intermediate 
regime arises. Several normalization conditions fixing the undefined constants are proposed, and estimations 
for the mass shifts are given in these cases. 
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I. INTRODUCTION 



There is a belief that quantum effects are only relevant for the black hole physics when the 
Compton wave-length h/Mc of the black hole is of the order of its gravitational radius IGMjc^ . 
This opinion is mainly based on the study of the Hawking radiation , which is incredibly small for 
regular black holes, or just on dimensional reasons. At the same time, we know that the metric field 
has a singularity on the horizon of a black hole, at least, in the Schwarzschild coordinates. So, we 
can expect a new physics in its neighborhood, which is not catched by dimensional considerations 
since there are large dimensionless parameters - the metric components. The aim of this paper is 
to derive the effective potential of the Higgs field on the Schwarzschild background, which governs 
the masses of particles in the standard model. It is evident that the effective potential must depend 
on a point of the space. Moreover, it is almost evident that it should be singular at the horizon 
of a black hole as the metric components are singular. So, we expect that the masses of particles 
considerably change near a black hole. 

In the seminal paper by Coleman and Weinberg H] , it was suggested that quantum fluctuations 
may change the effective action in such a way that a spontaneous breaking of the gauge symmetry 
occurs. That is, quantum fluctuations change a landscape of the potential entering the effective 
Lagrangian, and the properly constructed perturbation theory for a classically massless gauge 
invariant model will contain massive vector bosons. Investigating below the effective potential of the 
Higgs fleld on the Schwarzschild background, and, more generally, on an arbitrary static, spherically 
symmetric space-time, we shall see that the contributions to it due to quantum fluctuation are of 
the same form as in the Coleman- Weinberg potential, but with coefficients depending on a point of 
space. It turns out that, in interpolating the one- loop results, quantum fluctuations are so strong 
near the horizon of a black hole that they can restore spontaneously broken symmetry of the 
standard model or make massive particles to be very massive with inflnite masses on the horizon. 
More rigorously, if a conformal symmetry holds valid for massless particles on the considered 
background then there only exist two scenarios depending on a sign of the undetermined constant 
coming from a regularization procedure. In the flrst case, the symmetry is restored and all particles 
become massless at a flnite distance from the horizon of a black hole, while in the second case the 
masses of all massive particles grow, when one approaches the horizon, and tend to inflnity. These 
two scenarios can be called as the "hot" and "cold", respectively, as the matter falling into it 
becomes hot in the flrst case and cold in the second one. If the conformal symmetry is spoiled, 
the additional constant governing the mass shifts arises. Properly chosen, it can compel to tend 



2 



the vacuum expectation value of the Higgs field to a finite nonzero value on the horizon. So, 
the intermediate regime arises. It is clear that the mass shift results in many unusual physical 
effects and we do not enlarge on them. The effective potential also allows us to find how the 
energy density of the Higgs field varies from the point to point in a static, spherically symmetric 
space-time. However, to this end we have to fix additionally one arbitrary constant. 

Deriving the effective potential we, of course, shall make certain approximations. First, we use a 
toy model instead of the standard one. We consider a scalar SO{N) multiplet interacting with the 
gauge fields taking values in the Lie algebra so{N). The gauge symmetry is spontaneously broken 
by the Higgs mechanism down to SO{N — I). In the upshot, we choose = 4 in order to have three 
massive vector bosons as in the electroweak sector of the standard model. Second, we shall find 
the one- loop correction to the effective action only. In other words, we assume that the model is in 
a perturbative regime, i.e., the energies are sufficiently high and the coupling constants are small. 
This enables us to cast out the higher order terms of the loop expansion. Third, calculating the 
determinant emerging from an integration of the Gaussian path-integral, we exploit a quasiclassical 
(short-wave) approximation. That is, we take effectively into account the modes with the wave- 
lengths much smaller than the gravitational radius rg . The contribution of the long- wave modes to 
the effective action is well approximated by the renormalized classical action and we also add it to 
the final result. In fact, we add the kinetic term only neglecting a self- interaction of the Higgs field 
by means of the gauge fields. Notice that the quasiclassical approximation allows us not to make 
any distinction between the one-loop contributions to the effective action of one scalar particle 
and one polarization of a vector boson regardless of the fact that they interact differently with the 
gravitational field. Calculating the functional determinant, we do not use the standard method of 
the heat kernel expansion 0, ^. It gives an asymptotic series in the powers of 1/m, where m is 
a mass of the particle, which is proportional to the Higgs field for the vector bosons. The higher 
coefficients of this series are unknown, while to take a few first terms from it is unacceptable for 
our goal to obtain a regular effective potential of the Higgs field. 

The paper is organized as follows. In Sec. we give general formulae regarding the effective 
potential. We start with the BRST-quantization of our toy model and find the mass spectrum 
and free propagators in the Feynman gauge. Then, we obtain a formal expression for the one-loop 
effective action and reduce the problem to a finding of an analytic expression for Spln(V^ + m?) 



on the Schwarzschild background. Sec. HI is devoted to this problem, where it is solved in the 
quasiclassical approximation. Of course, we need some regularization prescription there to remove 
divergencies. This is done by introducing a certain cutoff parameter. It turns out that the obtained 
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expression for the one-loop correction is immediately generalized to an arbitrary static, spherically 
symmetric background metric. Also, in this section we investigate how the result depends on 
the regularization scheme applied. The form of a mass dependent part of the effective potential 
proves to be independent of the regularization scheme. Moreover, under the assumption that the 
conformal symmetry takes place in the massless case the form of the whole effective potential 
is fixed unambiguously. In Sec. we describe the renormalization procedure and provide a 
brief analysis of the physical implications following from the obtained one-loop effective potential. 
Possible scenarios of the physics near a black hole are established. In Appendix, we give asymptotic 



expansions of the integrals encountered in calculating the first quantum correction in Sec. ID. 

We use the following notation and conventions. The system of units is chosen so that the velocity 
of light c, the Planck constant h, and the Schwarzschild radius of a black hole are unities. In Sec. 



[V, where we analyze the physical effects near a black hole, we restore the ordinary length units. 
Greek letters denote space-time indices. Sometimes we shall use boldface characters to denote the 
spatial part of coordinates. Also, we use the mostly minus signature of the space-time metric g^i,. 

II. EFFECTIVE POTENTIAL: GENERAL FORMULAE 

Consider the model of the SO{N) scalar multiplet (j)'^{x) in the fundamental representation. 



n = 1,A^, minimally coupled to the gauge fields A^{x) = A^tJJ^, a = l,N{N — l)/2, on the 
Schwarzschild background g^u{x). The action functional of the model reads 



^{D,rf-Vo{cP)-\G%G^'' 



(1) 



S[(t>,A^] = j d''x^\g 

where Vo{(^) = A + ^l'^(|?/2 + \{(tPf/^, 

D^<t> = d^<p + ieA^cP, G% = d[^Al^ + ef^^A^Al, (2) 

and A, fi, A, and e are the coupling constants, /^^ are the structure constants of the Lie algebra 
so{N) with generators t^: 

[^aj^fc] = ^fah^C) = ta, tg. = —ta, (3) 

where transposition is understood with respect to the Euclidean metric preserved by the action 
of the SO{N) group. The Latin indices a,b,c. . . are risen and lowered by the Killing metric and 
the Greek indices /x, z/ . . . are risen and lowered by the space-time metric. Henceforth the matrix 
notation is implied and the indices are suppressed when it does not lead to confusion. The SO{N) 
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symmetry of the model is spontaneously broken by the Higgs mechanism down to SO{N — 1). The 
Higgs field can be decomposed as 



(x) = r]{x) +xix), 



(4) 



where the vector r]{x) is its vacuum expectation value. It provides a minimum to the classical 
potential Vo(</>) in the tree approximation. 

As the model (||) possesses the gauge symmetries, we need to use the BRST-quantization pro- 
cedure (see, e.g., |§, 0, §) in order to pass to the quantum theory consistently. The quantization 
procedure is quite standard, but it is usually given for a constant background. Therefore, we trace 
some its basic steps for our case. 

It is convenient to change the coordinates so that g^^ = 1 and g^^ = 0, and return to the original 
ones in the total action with ghosts. Introducing the momenta 



{</.„(t,x),7r"(t,y)} 



ni /^Oi 
a — ~"-^a 1 



C(x-y) 



7r„ = 



{Atit,^),Uiit,y)} 



(5) 



151 Vlfl 
where curly brackets denotes the Poisson brackets, we arrive at the Hamiltonian action 

1. 
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+ ^S(V,ni - ie7rta(t>) - + ^ A</'^V - Vo{(P) , (6) 



where the overdots denote the derivatives with respect to time and 



(7) 



As usual for the relativistic invariant models, Aq are the Lagrange multipliers to the first class 
constraints 

,'^(x-y) 



Ta := - icTTtacP, {T, (t , x) , Tf, (t, y)} = efMt,^)- 



Therefore, the BRST-charge is 

Q, = I dxy^\g\ 



(8) 



(9) 



where we introduce the scalar fermionic ghost fields (c". Pa) and (c", Pa) with the Poisson brackets 



{c'^it, x), Pfe(t, y)} = {e{t, x), n(t, y)} 



^?(x-y) 



ghc = ghc = — ghP = — ghP = 1. 



(10) 



Other fields are graded to have zero ghost number. We choose the gauge fixing fermion ip so it 
results in the Feynman gauge 



-1. 



(11) 



Adding — {fi, tp} to the Hamiltonian action (^) taken on the constraint surface, we obtain the total 
Hamiltonian action with the Lagrangian 



L 



H 

tot 



.\'^ + \DiC^D'<l)-V^{<l))+&^{Pa-Pa)+d'Pa-d'c^diPa + eC^ 
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(12) 

Substituting this action into the functional integral and integrating over the fields c, P, U and tt, 
we have (cf. M) 



L 



tot 



1 



1 



- ^o(?? + X) + -^xy^(j])X + Q^xdi,^ + - -(dy.T]A^r] - riA^d^ri) - ied^xA^X - e\A^^A^T] 

- i^xAlA^x + c\t,A^^d^ - xtlhm - Id^.Al^f^A'^A^'^ - ^l{f-^A'^A^'^f] , (13) 

where Vq := 3^ Vo/5(/)" and we take into account that the space-time is Ricci flat: R^,^ = 0. By 
the same reason we did not add the nonminimal renormalizable coupling ^Rcj)'^ to the action (|l|). 
Such a term is exploited in some inflationary models Even if the gravitational field does not 
satisfy the vacuum Einstein equations these terms will be irrelevant for our subsequent analysis. We 
shall make all calculations in the quasiclassical approximation (see below) and these contributions 
can be neglected. Upon integration over momenta a local contribution to the path-integral measure 
appears (see, for example, [|lO|), but as we shall see such contributions are irrelevant for our 
calculations. 

The first line in (|l^) contains the terms quadratic in the fields and it is only needed to find the 
first quantum correction to the effective action. It is not difficult to deduce the mass spectrum of 
particles. Bearing in mind that for our model 



(14) 



we conclude that there are 2(A^ — 1) ghosts, 4(A^ — 1) vector bosons, and (A^ — 1) Goldstone bosons 
with the squared mass 



m 



2 2 

e f] , 



(15) 



the one Higgs boson with the squared mass 

mjj = fi'^ + 3Xri^ , (16) 

and the massless fields: (A^ — 1)(A^ — 2) ghosts and 2{N — 1){N — 2) vector bosons. Under the 
number of particles we mean the number of their states. 

The value of the effective action at <j) = rj when the other fields are integrated out is given by 
general formula (see, e.g., [§) 

ir[r]] = j DxDADcDPe'^'°'^'^^^''^^'''^\ (17) 

where the integral assumes a summation over the vacuum, one-particle irreducible, connected 
diagrams only. In view of the last remark, the linear in the fields x terms entering ([l3|) 

can be omitted. The one-loop contribution to the effective action is formally given by 

zTiM = Spln(V2 + m^) - Spln(V2^ + m^^^;) - i Spln(V2 + m^j) 

+ - - ^) Spin - - - Spin V^^ + . . . (18) 

where the dots denote higher terms of the loop expansion, which are small as long as the renor- 
malized coupling constants e and A are small. Also we omit the contributions proportional to the 



trace of the identity operator. The Feynman propagators appearing in the logarithms (18) are 
specified in the standard way by the ze-prescription tit? — ie. The operators and V^'^ 

are defined as follows 

(19) 

v'f = vH'i + /^[5A(r;^j + 2r;:,aA + rf^^c,^] = v^^^ + g^^mv^) + 2r;^,aA], 

where the square brackets enclosing indices mean an antisymmetrization without one half. In the 
last equality, we used the fact that the Schwarzschild metric is the solution of the vacuum Einstein 
equations. Formula ([l^ ) is valid up to contributions from the path-integral measure. It can be 
shown Q that the path-integral measure respecting unitarity of the S'-matrix gives a one-loop 
correction to the effective action of the same form as ([Tsl), but with an opposite sign and the 
Feynman propagators replaced by the advanced Green functions. A general scheme how to define 
and calculate the S'-matrix on a curved background can be found, for example, in [11 1. 



In the next section, we shall calculate the traces entering (|^) and the analogous expression 
coming from the measure in the quasiclassical (short-wave) approximation. In this approximation. 
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the derivatives in (19) acting on the right are of the order of unity and the connections are of the 



order of h (see for details [|12|, Hence 

Vli" = V^6^ + 0{h). (20) 

This approximation is adequate for the modes with the wave-lengths much smaller than the gravi- 
tational radius r^. It can be interpreted as we effectively integrate over the fast modes only in the 
functional integral (p!7|). The correction to the effective action (18) due to the slow modes with 



the characteristic wave-lengths larger than rg is well approximated by the total action (13) with 
renormalized fields and coupling constants. The ghost fields can be excluded since they do not 
influence the dynamics of the mean fields x ^/i- Thus, the one-loop correction to the effective 
action in the quasiclassical limit reads 

•rn~ ^(M i (^' + ^')f 1. , (V^ + m^)^ (iV-l)(iV-2) (V^)^ 
i^iim ~ -- iV-l)Spln --Spin Spin— (21) 

where the subscripts indicate what the Green function should be taken. The main problem facing 
us is to derive analytical expressions for Spln(V2 + m2)^ + . 



III. ONE-LOOP CORRECTION 



In this section, we find the one-loop correction to the effective potential given by quantum 
fluctuations of one scalar bosonic field in the quasiclassical approximation on the Schwarzschild 
background. It turns out, the expression we shall obtain is valid for an arbitrary static, spherically 
symmetric background. 

To this end, we write 

-^Spln(v2 + m2)^-i J d^xy^\Y,\M^)\^H-e{0+m\x)), J d^xy^\\i;^ix)\^ = I, 

^ (22) 
where il}^{x) are the eigenfunctions of the wave operator and —£{$,) is its spectrum, ^ specifies 
the quantum numbers. The approximate equality appears since we assume that m? depends on 
the space-time point and so Eq. ( p2D is valid in the leading order in h only. The correction to 
this equality is proportional to a commutator of the wave operator with m?{x) and we neglect it. 
Inasmuch as we work in the Schwarzschild coordinates, the interval is 

ds'^ = (1 - r~^)dt'^ - (1 - r-^yUr^ - r^{de^ + sin^ 9d^'^). (23) 
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Besides, we assume for a while that ip^{x) vanishes at the big spatial sphere surrounding the black 
hole and is T-periodic with respect to the time variable. Hence, we have to solve the equation 



(1 _ r-Y^d^ - \dr{r^ - r)dr - . ^ ^ M nOde sin^a^ + dl) 



^ = -s^P, (24) 



with the given boundary conditions. Separating variables 

i^nkim = e'^%m{0, ^)Rnki{r), E = Ink/T, (25) 

where Yifn{9,(p) are the spherical harmonic functions, we reduce Eq. (|2^ ) to the radial one 

-^dr{r^ - r)dr - (1 - r'Y^E^ + ^^^-^il R^^l = SnklRnkl. (26) 

It can be solved in terms of the Heun functions (see, for example, |jl^), but the properties of 
these functions are not well investigated. The spectral problem ( |2^ ) was studied in many respects 
and the approximate solutions to it were found under certain assumptions (see, e.g., |Q, |l5|, [l6||). 
However, keeping in mind that we are interested not in the spectrum itself, but in (^), we shall 
try for another strategy and seek for a formal quasiclassical solution of ( [2^ ) at arbitrary values 
of the energy and angular momentum. We shall not write the Planck constant explicitly as it is 
restored in an obvious manner. 

Making the standard substitution, 

i?(r)=/(r)e^^W, (27) 
in (p6[) we arrive in the two leading orders in h at 



S' = ±(1 - wrWE' -{l-w){E + L^w^), f= (28) 

where w := r~ '\ := + and c is some constant. The turning points are at the the vanishing 
radicand 

E^-{1- w){enki + lW) = 0, (29) 

or at the boundary sphere or at the origin of coordinates. The spectrum is determined by the 
requirement that (^) is a single-valued function 

J dr{l - w)-^ ^E^ - (1 - w){£nki + L'^w'^) = vrn + 7, n G Z, (30) 

where the integral is taken between two turning points and 7 depends on whether the particle 
strikes the boundary sphere or reaches the origin or not. The values 7r/4, tt/2, and 37r/4 of 7 



correspond to the cases: (i) the particle moves between the origin and the turning point (p9|), (ii) 
moves between two turning points (M) or between the origin and the boundary sphere, and (iii) 



moves between the turning point ( p9D and the boundary sphere. It can be proved |12, |T3| that 
the spectrum derived in this way differs from the exact one by the terms of the order of h?. The 
normahzation condition ( p^ ) becomes 

T|c|2 r = 1, (31) 

where r\ and ri are the turning points and we use the standard normahzation of the spherical 
harmonics: 

y d^de sin QYtJe, ^)Yv^, {e, if) = 5u,5mm' . (32) 
Recall that the spherical harmonic functions possesses the property 

\Yue,v)? = '^-^- (33) 

m=—l 



\cnki?\n{-enki + rn^) 



Then the integrand of (^) can be written as 

-. ^ oo oo 

(34) 

In the quasiclassical approximation and in the limit T ^ oo whereas the radius of the boundary 
sphere tends to infinity, it is admissible to pass from a summation over the quantum numbers to an 
integration over them. The Bohr-Sommerfeld quantization condition ( pO| ) and the normalization 
l|) entails 

dn = -- r = (35) 



The Jacobian is nonzero everywhere. Thus, we obtain for the one- loop correction ( |34[ ) 

1 dE dL^ de e(E^ - (I - w)(e + L^w^)) ^ , 2s 
/= / / ^/ \ ^ ^ rln(-g + m^ , (36) 

where 9{x) is the Heaviside step function. Tracing all the steps of our derivation starting from (p^, 
it is not difficult to see that for an arbitrary static, spherically symmetric metric in the adapted 
coordinates we arrive at 



' ' ' ln(— e + m ), (37) 



Svr 2tt Jq 27r ^E^ - goo{e + L^w^) 
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where we assume that goo9rr < 0, and ^'s denote the Kilhng vectors 

6 := dt, ei := 9^.m'. + eyCy + e.C) = -2r^ (38) 

with being the generators of rotations. In passing, we cast the integral I into an exphcit generally 
covariant form. According to the ie-prescription for the Feynman propagator, the logarithm should 
be written as 

ln(-e + m^)F = ln| - e + - i'K0{£ - m^). (39) 

As for the advanced Green function we specify it by adding ie to the energy E. This is a correct 
definition of the advanced Green function, at least, for the space-time points above the horizon 
of a black hole. It is that region which we shall be interested in. Hence, for the advanced Green 
function the logarithm should be replaced by 

ln(-e + - ieg^^E) = ln| - e + - m s\gii{goQE)0{£ - w?), (40) 

in the quasiclassical approximation. 

The obtained integral (^^ has ultraviolet divergencies. Since this is a multiple integral it should 
be regularized with some care. In contrast to the case of a flat space-time, there are not continuous 
symmetries like the Lorentz symmetry which we should preserve. The only residual symmetry 
is the discrete T-symmetry E —E. Therefore, we regularize the integral cutting the range 
of integration over the respective variables preserving the time reversion symmetry. At that, we 
apply this regularization prescription to the initial integral (^) or more rigorously to the initial 
sum ( |3^ and assume that the regularization parameters do not depend on and r. In other 
words, we can find the derivatives of the effective Lagrangian ( |3^ ) with respect to and r by 
a mere differentiation of the expression under sums^. Then, making changes of variables, which 
do not depend on and r, we reduce this integral to the integral which is divergent in one 
variable only. That is, the "singular" boundary of the integration domain is retracted to a point. 
The regularization procedure of singular integrals in one variable is well developed (see, e.g., |p!7|| ). 
We simply cut off the divergence of this integral. The emergent regularization parameter is some 
function of the initial regularization parameters which does not depend on m? and r. An exact 
form of this function is irrelevant for us. All we need to know about it is that it tends to the 

^ Such a trick is a standard tool of the regularization on a flat space-time. Unfortunately, it is unapplicable in our 
case since the integral remains divergent after differentiation and we lose the information about the terms of the 
effective potential at lower powers of m^. They are not constant as in the flat space-time. 
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regularization removal limit as the initial regularization parameters do. Notice that the change of 
variable (|35| ) is independent of and r. 

Consider the contribution ( |37| ) to the effective Lagrangian at the point above the horizon of the 
black hole {w < 1). The integral over the squared angular momentum is easily taken with the 
result 



ln(— e + m^)F-, 



(41) 



where y = goo = ^f. Despite the real part of this integral does not contribute to the effective 
Lagrangian, because it is canceled out by the same contribution from the advanced Green function, 
we investigate its asymptotic expansion as well. In order to take the obtained integral, it is useful 
to pass into the polar system of coordinates 



£ = pcosip, E = psinip. 



(42) 



Then the integral over the angular variable converges, whereas the integral over p diverges 
quadratically. We regularize this divergence introducing the cutoff parameter A 

d(p ^1 — y cot if 



I 



167r3 



'arctany -lu/i y Jq 

After a suitable change of variables, we come to 



dpp[ln I — p cos if + m \ — iiT9{p cos ip — m )]. (43) 



Re/ 



1/2 



327r3y (1 - 2;2)3/4 



A2 

In A + A^ln 

2 



m 



In 



m 
xA 



+ 



xA 



and 



Im/ 



m 



dx 



yjl-xy /A 



327r2y 7(A2/m4_i)-i/2 1 + 



x^ 



(44) 



(45) 



To obtain an asymptotic expansion of the real part of (44), in terms of the regularization 
parameter we use the formulae (A5) and ( |A11 ) given in Appendix. Then its singular and finite 
parts look like 



1 rA2, a2 

e 



Rel= In— /o(y) + A2/o(y)-2Am2 /i(y) - In ^1 + 



m 



h{y) + 2h{y) - 3 - 3^1 + - | ^/l + y^ - In ^/l + y^ - + ^1 + ^2 in(l + y^) 



+ 21n^ (l + + - + - h{y) 



., (46) 
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FIG. 1: At the left: the function goiy) at /3 = l/ln2 (thick Une), at /3 = 1.4 (dashed Une), and at /3 = 1.34. At the 



right: the function g2(y)/7r (thick hne), the function gi{y) (dashed line), and ~(R{y)/ny). Bars over the functions 
mean that we subtract from them their values a.t y — 1. 



where 
Io{y) 

io{y) 
h{y) 
h{y) 
h{y) 



dxZ{x) 



dxZ{x) In |x| 



Try 



dx- 



dx 
dx 



X 

Z2{X) 



y 



2 ./_oo 

dx 

-oo X{x'^ 

y 



In- 



nf-Z yZ + 



x/(l - + y^) -y 



1 + Vi + y'-^in^i + y', 



^1^'^ In Ixl 



dx 

2 ./— oo •^^ 



+ 



Vx2 + y2 ' 2 x2 + y2 
All the integrals on the left hand side are taken over x G [—1, (1 + y'^)"^^'^] and 



In 



x2 + y2 



Z(x) 



(^/r^^-xy)i/2 



1 - — + 0(x ). 



(47) 



(48) 



(1 - x2)3/4 

The integrals Iq and I2 are reduced to dilogarithms, while Ii is an elliptic integral. We see that 
the term at In m? in Eq. ( ^6|) vanishes. Introducing new notations for the infinite constants 



entering the real part (46), we can write 



Re I = a 



m 



0- 



y 



+ Bi — 

y 



/i(y)-lnVl + y2 



m 



32n^y 



i2(y) - 1 + VT+^ (in yrr^ - 1) - i (i"' ^1+^ - 



= : Bogo{y) + Bivn^gi{y) 



m 

327r3 



92 (y), (49) 



where /3 is some dimensionless constant. 
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The integral (^) is easily expanded in (A^/m^ — 1) So, up to the terms vanishing at 

A ^ oo, the one-loop contribution to the effective Lagrangian reads 

, R(v) . , ^ / 9s , X 

= ^o^^ + ^i— +^2^4-— ^In ym^ , 50 
y y D47r^ 

where 



- Y^v'l + 2/2-^arccothyiVl + y2 + ^- (51) 

The analogous contribution ( pO|) from the advanced Green function is zero. Recall that the obtained 
expression is valid for the points above the horizon of any static, spherically symmetric metric. In 
spite of we use the system of units in which a mass is measured in hr~^c~^, formula ( |50| ) holds as 
is in the system of units H = c = 1 with the standard mass units (say, GeV/c^). The undetermined 
dimensional constants are only redefined and y = 1 — Vg/r. The plots of the functions gi{y) and 
R{y)/y are presented at Fig. |l[ Asymptotics of these functions are as follows: at the horizon 

R/y = -^v'^ + 2 1^'' I " ^) + 16^ + 

50 = ^(1 - /31n2)y-i + ^(1 - /31n2 + + o(y), 

8 (52) 

ffi =lny-l-31n2 + o(l), 
92 = -vry"^ + 0(2/"^), 
while at the spatial infinity 

R/y = 0.561 - 0.971(y - 1) + 1.23(y - 1)^ + 0((y - 1)=^), 

50 = 3.45 - 2.10/3 - (2.95 - 2.19/3)(y - 1) + (3.04 - 2.22/3)(y - 1)^ + 0{{y - if), 

(53) 

51 = -3.04 + 1.02(y - 1) - 0.534(y - 1)^ + 0{{y - if), 
92 = -0.313 + 3.45(y - 1) - 3.20(y - 1)^ + 0((y - 1)3). 

We provide the approximate numerical values for the expansion coefficients at the spatial infinity. 
Their exact expressions through the integrals or arccotangents of huge arguments can be read off 
from Eqs. (|4^), 

Let us study how the obtained asymptotics depend on a shape of the integration domain in 
the space of the quantum numbers. The asymptotic expansion of the imaginary part (Ea) for an 
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arbitrary star-shaped region of integration with the boundary A(x), A(2;) oo, takes the form 



reff ^ r dx yi-xj/ 2. s A(0) 
1 ' 327r2y 1 + x2 16^2 y 



m 



4 



327r2 



lf|„„|4A(0)lV.-^'<°> 



4 



647r2 



In(ym^), (54) 



2 V2 ' ' "J A(0) 
where x = corresponds to massless modes. So, the part of the one-loop contribution to the 
effective Lagrangian depending on a mass is completely determined by three constants irrespective 
of an integration domain, which, of course, should coincide with the whole space of quantum 
numbers in the regularization removal limit. The massless contribution depends considerably on 
the shape of the integration region. We can fix the form of A(x) unambiguously if demand that 
in the massless case the effective action possesses a conformal invariance on the given background. 
From Eq. ( |37|) and ( |4l| ) it is not hard to see that a stretching of the metric by a function A(r) is 
equivalent to a stretching of y by this function. Therefore, if we take A2(x) = (1 + a;2)A2(0), the 
one-loop contribution to the effective action, 

4 

^ = Aoy~'^ + Aim^y-^ + A2m^ + A2m^y-^ - In(ym^), (55) 

will be conformal invariant in the massless case. The factor A^2 canceled by the factor coming 
from the measure. The choice A2(x) = (1 -|-x2)A2(0) corresponds to a usual energy cutoff, i.e., all 
the modes with the energies higher than A^^'^{0) are thrown away. In this case, the additional, as 
compared to ([50|), constant A2 arising in Eq. ( |55[ ) must be zero. 

As far as the real part of ( ^7|) is concerned its dependence on a mass, and, in particular, 
vanishing of the term at m^lnm, is not influenced by a changing of a shape of the integration 
domain. However, the functions gi{y) change with A(x) analogously to the first term in ([5^). 

In the next section devoted to renormalization we shall analyze the one-loop contribution to the 
effective action in the form (^5|) keeping in mind that a spoiling of conformal invariance modifies 
the massless contribution. 

IV. RENORMALIZATION AND ANALYSIS 

In this section, we describe the renormalization procedure and give a brief analysis of possible 
physical effects following from the form of the obtained effective potential. A thorough study of 
the physical implications requires a detailed investigation of experimental data, modification of the 
obtained results to the standard model instead of the toy model (|l]) we used, and comparison with 



other approaches (see, e.g., [||, 16, |T^]). It will be given elsewhere. 
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First of all, we take into account that the effective potential should reduce to the Coleman- 
Weinberg potential on the spatial infinity y = 1. This, in essence, results in a subtraction of 
values of the functions of y in the effective Lagrangian (^) taken at y = 1, and replacing all the 
coupling constants by their renormalized values. We shall not mark them by new labels and just 
bear in mind that such a renormalization has been done. Notice that the undefined constant A2 is 
completely absorbed by this renormalization. The residuary constants Ai have to be determined 
by some normalization conditions and, eventually, should be taken from experiments. As we have 
no any data on their values, we shall study possible scenarios for different Ai. 

We start our investigation with the case A2 = 0. Then, on substituting the masses of the 
particles in our model (^), (p^) expressed in terms of tj, the one-loop effective potential reads 

V = A + N^Aoil - + - y-i) + ln(y|/.2 + 3Xrj^\) 



2 J 



l + ^ln(y|/.2 + 3A7?2|) 



Ai[3(iV-l)e2+3A](l-y-i)| 



,(X 3(iV - 1)6^ + 9A2 3(A^-l)eMn(eV) + 9A2ln|//2 + 3Ar;2n 
+ + J' ^''^ 

where A is the cosmological constant minus the value of the residual part of the potential taken 
in its minimum at the spatial infinity. For simplicity, we assume that the cosmological constant 
vanishes. Besides, the logarithms of the energy scale at which the coupling constants and masses 
are measured are included into A, ^u^, and A. Because of that the expression for the effective 
potential ( |56| ) and some expressions below look wrong by dimensional reasons, but it is just an 
appearance. 

In order to make our toy model being similar to the electroweak sector of the standard model, 
we choose = 4 and take the following values of the coupling constants and masses [^] at y = 1 

2 

mo = 83, ruHo = 150, % = 247, = -12445, e = — , A = (57) 

Vo 2r?^ 

where the value of fj? is taken so that the effective potential at the spatial infinity has a minimum 
at the experimental expectation value of the Higgs field tjq. Of course, choosing A^ = 4 we have 
too many massless particles as compared to the electroweak sector of the standard model, but the 
difference in their number is irrelevant for our qualitative analysis. 

On the Schwarzschild background, the one-loop effective action has an extremum when the field 
ri(y) satisfies an analog of the Ginzburg-Landau equation 

2 



(58) 
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where we already regard rj as the one component of the vector rj"', other components being zero. 
Besides, we neglect a self- interaction of the field r] through the gauge fields ^4^. The kinetic term 
arising in ( [5^ ) takes effectively into account the long-wave modes. Below we shall see that it is small 
on the functions minimizing the effective potential V{7]) in a sufficiently wide range of changing 
of the variable y and can be considered as a perturbation. In that case, the below analysis is 
immediately applicable to an arbitrary static, spherically symmetric gravitational field. 

In the conformal case, A2 = 0, a crucial role for the physics near the black hole plays a sign of 
the constant Ai. Indeed, ignoring the kinetic term, we have 

2r] 



V 



dyV 
'V"{rj) 



yV'iv) 



y 



[3{N - l)e' + 3A] + 



3(Ar - l)e^ + 9A^ 
32^ 



3A + (iV - l)e4A-i 



(59) 



The second term in the curly brackets is positive at the spatial infinity, because it vanishes at so 



small r] that the Higgs boson has an imaginary mass according to (16). Consequently, if > 0, 
the function rjiy) is monotonically decreasing since V" > at the minimum. 

Assuming r] is finite at y ^ 0, it is not difficult to obtain from Eq. (|58| ) without the kinetic 
term that either r] tends to zero or 



— [3(iV-l)e2+3A] 

y 



3(iV - l)e^ + 9A2 
32^^2 



1^ 



3A + (iV - l)e4A- 



Iny. 



(60) 



Therefore, Ai should vanish, so does the term on the right hand side, but it is impossible as we 
have already seen above. Hence, the only two possibilities exist: r] tends to infinity or to zero at 
the horizon of the black hole. 

At Ai > 0, the function r]{y) is monotonically decreasing and tends to infinity at y — > 0. To 
understand a behavior of r/ at the negative Ai, we rewrite the minimum effective potential condition 
under the assumption y ^ and 77 ^ 00 in the form 



^i[3(iV-l)e2 + 3A] 



+ yv 



A 3(A^- l)eMne^ + 9A^ln(3A) 3(iV - l)e^ + 9A^ / 1\ 



0. (61) 



2 ■ 327r2 ■ 327r^ - / ■ 2 

Hence, r] tends to infinity as y^^^^. It is not hard to investigate the obtained asymptotic equation 
on yrj'^. At Ai < 0, it has two solutions, which coincide at Ai, only if 

1 {N - l)e^ + 3A2 



Ai > Ai 



327r2 (iV-l)e2 + A 

3 IQn^X + 3{N - l)e^ In + dX"^ ln(3A) 
'2 3(A^ - l)e4 + 9A2 



X exp 



-4.66 • 10"^-^. (62) 
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The right dimension can be recovered if we recollect that A contains the logarithm of the energy 
scale at which the coupling constants are measured in experiments and take the values (|57|). The 
value of Ai depends on this scale through the coupling constants. Thus, at Ai G [^i,0], the Higgs 
field goes to infinity at the horizon, while at Ai < Ai it vanishes there. 

Assuming r] to be small, we can expand V'^rj) around the point i] = with the result 



V = i) |3(A.-l)e^ + 3Al + f [l + ^ (ln|,,^| + i) 



+ \ / r,^W + ... 63 
327r^ 



We see that at Ai < and the sufficiently small y there exists a value of the function r] that the 
written terms disappear. Equating the first two terms in (^) to zero, we can determine the critical 
point 7/0; where r] vanishes. The derivative of the Higgs field rj'd/) at this point is infinity. In the 
narrow interval Ai £ {Ai, 0), the real part of the effective potential has three extrema, two of them 
being minima. That is, in this interval the system may be found in a metastable state at small 
y. Notice that all the described features are typical for the Landau phase transition theory. As 
for the kinetic term, which we neglected in our considerations, it merely accelerates the process of 
going rj to zero or infinity. 

So, we distinguish two regimes which can be called as the "hot" and "cold". In the "hot" 
scenario, Ai < 0, and the spontaneously broken symmetry is restored near the black hole. All 
particles falling on it become massless, the matter is disrupted and electrically charged particles 
radiate a huge number of high-energy photons which freely from the boundary of this region reach 
a remote observer. In the "cold" regime, Ai > 0, the masses of all massive particles grow, when 
they approach the gravitating body. The high energy physics is suppressed by large masses, the 
matter rapidly gets cold and shrinks. The Bohr radius decreases, so does the range of action of 
the nuclear forces, while the energy level spacing increases. 

Speculating in this way, we interpolate the one-loop results to a region where quantum cor- 
rections become large. Furthermore, the large gradients of rj produce a considerable gauge fields 
A^ which act on the Higgs field and change it. Thus, we can believe in the obtained implications 
to a certain extent only, expecting that we catch a qualitative behavior of the Higgs field near a 
stationary, spherically symmetric gravitating body. 

Now we turn to the case of a spoiled conformal invariance for massless modes A2 7^ 0. It is 
not difficult to see even from ( [5^ ) that at A2 > the expectation value of the Higgs field grows 
exponentially near the horizon. Its asymptotics reads 

y 327r^ 
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This is a highly unhkely scenario. At ^2 < 0, we have two regimes: 

Ai 



2A2 

where 



(65) 

1' 



A{ := ^-^^ ^ 4.39 • 10^. (66) 

^ 3(iV-l)e2 + 3A ^ ' 



In the first case, the phase transition occurs at the critical value uq determined by the equation 

(67) 



6A ( 1 - - ) (Ai - 2A\A2) = A\ 

yo. 



3A / , 2, 1 



The expectation value of the Higgs field has an infinite derivative in this point with respect to y. 
In the second case, the Higgs field tends to the finite nonzero value on the horizon 

. _ 3(iV - 1)6^ + 3A / A, \ 
(0) - s{N - l)e^ + 9A2 - ■ 

For example, at vanishing A1/A2, quantum fluctuations of the Higgs field become tachyonic there. 
Thus, a breaking of the conformal invariance for massless modes may give rise to stabilization of 
the masses behaviour near the black hole. 

A. Normalization conditions 

As we have already mentioned, the constants Ai should be determined from experiments. How- 
ever, we can fix some of them demanding a fulfilment of certain properties from the Higgs field rj 
or the effective potential V{ri). As the effective potential V is regular in y at y = 1, the field r]{y) 
is regular as well. Hence, r] approaches a constant value at the spatial infinity as 1/r. It is possible 
to choose the constant Ai so that rj tends to a constant faster, i.e., it goes to rjQ as l/r^, keeping 
A2 = 0. From Eq. (^) we deduce the corresponding value 

Acr_ 3A/i2 + [3(jV-l)e^ + 9A2]r?g ^ 

~ 327r2[3(iV-l)e2 + 3A] ^^^^ 

That is, we get to the "hot" case. If yli < and Ai > Af then the masses initially grow, when 
one approaches from the spatial infinity to the black hole, up to certain maximum values and then 



decrease down to zero. If Ai < Af, they monotonically tend to zero. For the marginal value (69) 
the Higgs field behaves at the spatial infinity like 

77 « ryo - 0.562(2/ - 1)^ + . . . (70) 
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In particular, assuming that the gravitational field in the outer space of the Sun changes slowly 
and is spherically symmetric, at least approximately, the relative mass shift on the Sun's surface, 
where y — 1 ~ —4.24 • 10^^, becomes 

^ « -4.09 • 10-1^. (71) 



The critical point yo determined by Eq. ( |6^ ) or (p7\ ) is found to be 

yo « 9.30 • 10"^ (72) 

Recall that the kinetic term makes yo to be larger. For comparison we give the analogous quantities 
for the "cold" regime at Ai = 0: 

?7f«r/o-1.12(y-l) + 0.574(y-l)2 + ... ^ 1.92 • 10"^ (73) 

Vo 

Assuming A2 7^ 0, we can make the masses of particles to be almost constant. For example, 
demanding the first two derivatives of ri{y) to be zero at the spatial infinity, we find 

Ai « 2.89 • 10^°, A2 « -1.21 • 10^ ^^"^ ~ ^° « -4.80 • 10"^°. (74) 

■no 

However, the price we should pay for stability is the large energy density of the Higgs field, i.e., 
the large value of the effective potential taken in its minimum, which, moreover, is not integrable 
at the spatial infinity. 

Another reasonable normalization condition consists in the requirement that the energy density 
is integrable at the spatial infinity. To meet this requirement we should expand the potential in a 
Taylor series around the point y = 1 up to (y — 1)^ order and assign its corresponding coefficients 
to be zero. At that, we have to exploit Eq. (^) without the kinetic term to express the derivatives 
of the Higgs field. After a little algebra, we arrive at the three equations on the three unknowns 
^0) ^1; and A2. They have two solutions 

i) Ao~-9.36•10^ AlRi3.58•10^ ia 1.04 • 10-^ 

(75) 

ii) Ao ^ -1.01 • 10^ Ai a 4.24 • 10^, A2 « -3.80 • IQ-"^. 

In the first case, masses of all massive particles grow exponentially in approaching the horizon. In 
the second case, the expectation value of the Higgs field decreases and tends to the finite nonzero 
value 

^fc^ ^ -0.220. (76) 

Vo 
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FIG. 2: The dependence of a mass of the vector boson on the metric component goo = y- The smaU plots separately 
depict the corrections to the mass due to the kinetic term regarded as a perturbation. At the left on the top: the 
constant Ai ~ —37.9 as in Eq. (^) and A2 = 0, the dashed line corresponding to the mass dependence with the 
kinetic term omitted. The discontinuity is at the point of a zero Higgs boson mass, where the effective potential is 
not smooth. At the right on the top: the constants Ai and A2 are taken to be zero. At the left on the bottom: the 
constants Ai and A2 are given in the first line of Eq. ([tsI). At the right on the bottom: the constants Ai and A2 
are given in the second line of Eq. (JTq) . 



The asymptotics at the spatial infinity and the relative mass shifts on the Sun's surface are as 
follows 



It 



i) rj 

7] 



r/o - 107(y - 1) + 84.0(y - 1)2 + . , 
r]o + 13.3(y - 1) - 10.4(y 



1? + 



V 
V 



-1.4 -10^(2/ -1)4 + .. 
3.65 -10^(2/ -1)4 + .., 



^..1.84.10- 

Vo 



-2.28 • 10" 
(77) 

In spite of the requirement of integrability of the energy density at the spatial infinity entails 
that the conformal symmetry should be spoiled for massless modes, we assume that the massless 
contribution to the one-loop effective potential is well approximated by the conformal invariant 
expression, which we derived. If the conformal symmetry for massless modes strictly holds, the 
energy integral of the Higgs field diverges at the spatial infinity, at least, logarithmically. 

The plots of the masses m{y), which are proportional to i]{y), and corrections to them due to 
the kinetic term in all these three cases are presented at Fig. |2|. 
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APPENDIX A: ASYMPTOTICS OF SOME INTEGRALS 



In this appendix, we consider asymptotics of the integrals, which we encounter calculating the 
effective potential. We need an asymptotic expansion for the four integrals: 



.6 

Li = dxip{x)lii\x — v\, 
Jo 



L2 
Ml 
M2 



dxcp{x) ln(x + v), 



.6 

/ dx(p{x) 
Jo 



Inll - ^1 + ^ 



(Al) 



X 

dx(p{x^ ^ ~^ ^ 





X\ X 
V 



X^ 



where v tends to zero from the right and (p{x) is an integrable on [0, b] function smooth at the 
point X = 0. The first two integrals should be expanded up to v'^ and logarithmic corrections to 
it, while in the expansion of the integrals Mi and M2 we should retain the terms up to the finite 
part only. 

The first integral reads 

pv rb 

Li = dxip{x) ln(u — x) + I dxtp{x) ln(x — v) 
Jo Jv 

rl rb—v 

= v dx(p{vx)[lnv + l'n.{l — x)] + / dx(p{x + v)lrLX. (A2) 
Jo Jo 



Expanding this in v we arrive at 



/•b r pb 

Li = dx(p{x)lrLX — V (^(0)(1 — Inz;) + (/5(0) Ine + / dx 



V 

y 



dx TT- 



+ o{v^), (A3) 



where e tends to zero. Analogously, for the second integral we have 



/•h r />6 

L2 = dx(f{x) In X + V </?(0)(l — Inf) + (^(0) Ine + / dx 



(fix) 



V 

T 



1 



Inv 



+ (p'{0)lne + J dx 



ip{x) 



+ o{v^). (A4) 



2 / s 

The obtained parts of the asymptotic expansions of the integrals Li and L2 pass into each other 
by the substitution v ^ —v and leaving \nv unchanged. Therefore, 



nb pb r p' 

L := dx(p{x) ln\x — v\ = / dxip{x) In \x\ — v 

J a J a \_J a 



X 



V 



J a 



dx 



ip2{x) 



+ <^'(0) In 



-cpiO) 



b 
a 

1 1 
b a 



+ 0{v^), (A5) 
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where a < 0, 6 > 0, and 



fc-i 



L J nn\ 



n=0 



(A6) 



For the integral Mi we have 



Ml 



dx(p{x 







dxv^(t;x)^ + / dx[^2(rE) + (/3(0) + /(0)x] / (A7) 



The second integral can be written to the required accuracy as 



L 



{x + vY V 



dx 



^2{X) , X 



In - +o(l), 







^/^^ , ln(x - 1) + X ^(0) 
dx 7, = 



x^ V 



\ V V / ^ V 



b/v 



99' (0) / dx 



ln(x — 1) + X 



X 



b 1 o V 



TT 



+ 0(1), 



+ 0(1). 



(A8) 



The first integral in (A?) is easily expanded in v. Collecting all the terms together we obtain 



Ml 







dx-^-i^ 0) (1 + ln- 
X V 



+ 



h 

f2{x 



\ mix) 

dx — 
x^ 

o 

TT 



+ rA^ln2_+ / dx^-^lnx- — (/p'(0) - ^(l + ln&) +0(1). (A9) 

Q X O 



By the same way, 
1 



Mo 







^x^-V.(0)(l + ln^ 
X V 



+ 



99(0) 



^x^ 



+ 



In"^ - + / dx 



.'^2{x) 







TT 



X^ 



Inu 







lnx+ — (/9'(0) - -^(1 + ln6) + 0(1). (AlO) 
6 b 

The expansions of the integrals Mi and M2 are very similar but not the same. Adding them, we 
can write 



M := 



dxip{x) 



ln|l-^| + ^ 1 



V I V 



'dx^+min 



x 



Inv 



dx 



(P2{x) 



+ ip'iO) In 



^(0) 



1 1 

b a 



+ l\x^ln\x\ 



x^ 



+ 





b 








a 



ln|H-^2)-VP(0)(--- + — 



1 1 In 6 Inlal 



a 



+ 0(1). (All) 



Note that the terms proportional to In v are canceled out. 
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